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We study some quantum theories which are divergent but for which renor- 
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malization can be performed nonperturbatively and explicitly. The result is 
a well-defined, finite formulation of these theories, in which neither a cutoff 



nor a bare coupling constant appears. Such theories describe 'contact' inter- 
CD ! 

actions between particles which are encoded into the boundary conditions of 



the wavefunction rather than in the hamiltonian. It is the attempt to de- 



scribe them in terms of conventional potentials or self-interactions that lead 
to divergences. 

We discover that, after renormalization, the dynamics is described by 
a new operator (the 'principal operator' $(i?) ): the Schrodinger equation 
for energy levels is replaced by the eigenvalue problem Q(E)if) = for this 

operator. (More generally, the resolvent of the hamiltonian has an explicit 

1 Department of Physics and Astronomy, University of Rochester, Rochester, New York 
14627, USA 



formula in terms of the inverse of $(£').) Moreover the interactions are not 
described by the boundary conditions (domain) of the operator, but by the 
formula for it; i.e., by its action on smooth functions. Even when the theory 
cannot be explicitly solved, it can be given a finite formulation once this 
operator is determined. It is then possible to apply standard approximation 
methods; in particular we can determine the energy of bound states, which 
would be impossible in petrurbation theory. 

We propose to call such theories which are apparently divergent, but have 
a finite formulation in terms of the principal operator, transfinite quantum 
theories. We construct some examples of such transfinite quantum field the- 
ories: quantum mechanics with 'contact' interactions, three body problem 
with contact interactions, quantum fields (fermionic and bosonic) interact- 
ing with a point source, many body problems with contact interactions, and 
non-relativistic field theory with polynomial interactions . 

As an application we develop a theory of self-interacting Bose fields in two 
dimensions with an attractive self-interaction. The ground state is a Bose- 
condensate for which the conventional many body theory breaks down due to 
divergences. The magnitude of the ground state energy grows exponentially 
with the number of particles, rather than like a power law as for conventional 
many body systems. 
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Keywords: 



1 Introduction 

There is a fact about classical mechanics that is so basic that we dont even 
call it a law of mechanics. Perhaps it should be called the zeroth law of 
mechanics: 

All particles behave like free particles for short enough intervals of time. 
In the language of modern analysis, the path of any particle is differentiable, 
so that it can be approximated by a straight line for small enough intervals 
of time. There are other physical phenomena (such as diffusion) where this 
law does not hold. Hence it is not a self-evident fact, but rather a law of 
nature. Newton's second law states that the deviation of the path from that 
of a free particle is second order in time (when the path is viewed as a curve 
in configuration space) and is given by the force divided by the mass. 

In quantum mechanics we have a similar zeroth law as well. The propa- 
gator (i.e., the integral kernel h t (x, y) =< x\e~ Ht \y >, where H = p 2 + V(x) 
is the hamiltonian operator) is a gaussian upto power law corrections: 



h t (x,y) = —— r ll + 0(t)\. (1) 

(47rt)2 V / 

(Of course this statement, just like the previous statement in classical me- 
chanics, is true only for non-singular potentials. The quantum theory is 
more forgiving of singularities than the classical theroy. We will consider 
some singular situations later.) 

Now let us look at quantum field theory. The analogous statement would 
be that the correlation functions G(x,y) =< <p(x)(j)(y) > of a quantum field 
4> will approach those of a free field theory when the distance \x — y\ between 
points becomes small. This is however not true for realistic quantum field 



theories. In Quantum Electrodynamics for example, quite the opposite is 
true: the smaller the distance between points the larger the deviation of a 
correlation function from that of the free theory. Asymptotically free theories 
(such as Quantum ChromoDynamics) are better behaved: their correlators 
do approach those of the free theory but the deviations are of order lo } , 
rather than \x — y\. 

Thus if we apply conventional ideas of dynamics to such quantum field 
theories we will run into divergences. These have to be removed by the 
unwieldy procedure of renormalization. This leads us to ask if there is another 
way to formulate such theories: to describe the dynamics of the theory by 
some operator (instead of the hamiltonian) which is able to describe the 
logarithmic deviation from the free theory. We must first construct such 
a new picture by working within the traditional renormalization method 
applied to simple model problems. Once we learn the basic ideas we might 
be able to deal with more realistic situations. In this paper we will propose 
such a reformulation of renormalization by studying a series of examples, 
starting from quantum mechanics (with singular potentials) and ending with 
some non-relativistic field theories. 

Such a new point of view is necessary, since the most fundamental theory 
of nature-the standard model of elementaty particles- is a divergent quantum 
field theory. It was discovered in part by the requirement of renormalizabil- 
ity: that the divergences can be removed order by order in perturbation 
theory by redefining the parameters. The main remaining problems in the 
standard model are due to the unpleasantness left over from this procedure. 
For example, the 'naturalness' problem in the Higgs sector is an extreme 



sensitivity (quadratic dependence ) of observables (like the Higgs mass) on 
microscopic parameters. 

There are also more practical reasons to look for a new formulation of 
renormalization. One of the main obstacles to solving for the dynamics of 
non-abelian gauge theories is also that its quantum theory is divergent. We 
only know how to remove these divergences by perturbative renormalization. 
But perturbation theory cannot describe the formation of bound states: a 
serious problem since all the observable states of the theory are expected to 
be bound states (confinement). 

There are also motivations that come from outside of particle physics. We 
don't yet have a quantum theory of gravity due to the divergences that arise 
in quantizing general relativity: this time the infinities cannot be removed 
by renormalization. Perhaps a better understanding of renormalization will 
lead to a way of quantizing non-renormalizable theories as well: there might 
a be a non-trivial fixed point for the renormalization group. 

String theory is a serious candidate for a quantum theory of gravitation 
and, possibly even a unified theory of all forces. Moreover it is a finite 
theory. Thus it may be the ultimate solution for many of the problems noted 
above. However, it should be possible to account for the spectacular success 
of renormalizable quantum field theories at current energies without having 
to resort to a new theory of spacetime at Planck energies. There should be 
a consistent nonperturbative formulation of quantum field theories whether 
or not the ultimate theory is a finite string theory. To draw an analogy, it is 
possible to have a mathematically consistent formulation of continuum (fluid) 
mechanics even though the ultimate description of fluids is in terms of a finite 



(but enormously large ) number of particles. Quantum field theory should 
have a self-contained formulation whether or not it is just an approximation 
to a more fundamental theory. 

We also mention another great success of the idea of renormalization: the 
divergences that afflicted the theory of second order phase transitions were 
finally understood by Wilson's development of the renormalization group. 
New ideas on renormalization might also deepen our understanding of phase 
transitions. We will give some examples in this paper. 

Thus, at the heart of each of the fundamental problems of modern theo- 
retical physics is a divergent quantum theory. Progress is not possible with- 
out a better understanding of these infinities. The situation is analogous to 
knowing the basic laws of mechanics but without knowing calculus: anything 
beyond the simplest examples are inaccessible. Indeed it was a deeper under- 
standing of the nature of infinity (Cantor's work on the transfinite numbers) 
that made modern analysis possible and by extension the modern theory 
of dynamical systems. We must look for examples which will help us de- 
velop such a general theory of divergent quantum systems: a truly analytical 
quantum mechanics. 

We will start with simple quantum mechanical models with singular po- 
tentials and progress to systems with several degrees of freedom, the 
three-body problem and then to non-relativistic field theories. Eventualy we 
hope to formulate QCD in this way, but even the cases we have studied so 
far address some physically interesting problems. 

To be more precise, we have found a way to reformulate some quantum 
theories that are divergent in the usual formulation in a new way that is 



manifestly finite: there is no need for a regularization and all the parameters 
are physical, not 'bare' coupling constants. We call such systems transfinite 
quantum systems. In these systems, interactions are not specified by poten- 
tials or any other simple modification of the hamiltonian operator: indeed 
on smooth wavefunctions the hamiltonian acts just like the free hamiltonian. 
It is the attempt to shoehorn these system into a conventional description 
in terms of a hamiltonian such as H = — V 2 + V(x) that leads us to diver- 
gences. The interactions are encoded not into the formula for the hamilto- 
nian but into the boundary conditions on the wavefunctions (i.e., domain of 
self-adjointness of the hamiltonian). Thus the resolvent of the hamiltonian 
(which is a sort of Green's function) contains the complete specification of the 
system. The resolvent has the complete physical description of the system, 
in terms of physical and not bare parameters. 

Upto this point the ideas are not very new: the resolvent has been worked 
out in many simple cases, most notably by Krein and his school || . Unfortu- 
nately in most cases, determining the resolvent is the same as exactly solving 
the system. We need a way to think of such 'contact' interactions which is 
free from divergences and yet does not require us to solve the whole dynamics 
first: we just need to solve the 'singular part' of the dynamics. Our main 
discovery is that there is a new operator, which we call the Principal Op- 
erator <&(E), which describes the dynamics of transfinite quantum systems. 
It is free of divergences, can be determined explicitly (usually as an integral 
operator) and is quite simple in many cases. There are no subtleties in the 
definition of its domain: the interaction is described by a term in the formula 
for <&(E). In this way the formulation of transfinite quantum systems is no 



more difficult than that of finite quantum systems. 

The eigenvalues of energy are given by the solutions to $(E)\ip >= 0; 
the scattering amplitude is determined by the inverse of $(i?). Of course 
the solution of the eigenvalue equation or the inversion of $(i?) is a difficult 
dynamical problem: as difficult as solving the Schrodinger equations would 
be in finite quantum systems. But the point is that it is no more complicated 
than that: we can apply the standard methods such as variational principles 
or perturbation theory to the principal operator, since now we have a formu- 
lation free of divergences. The system need not be exactly solvable for it to 
have a transfinite formulation. Formulating the system amounts to finding 
the principal operator $(i?) while solving it amounts to inverting $(E"). 

The principal operator can be thought of as an effective hamiltonian 
obtained after 'integrating out' (or eliminating) the short distance degrees of 
freedom. This is reminiscent of Wilson's program of renormalization. The 
main new point is that we can get a closed experession for such an operator, 
rather than give an implicit prescription to find it as in Wilson's program. 

Simple, exactly solvable , examples of such transfinite quantum theories 
have been known for a long time: the two dimensional delta function po- 
tential is a good example ||. We will in fact start our discussion with such 
a simple example, the pole model of low energy scattering in quantum me- 
chanics. It has been known for a long time that at momenta small compared 
to the size of the scatterer, the scattering amplitude of a particle tends to 
a universal form f(k) = £ -i_ ik , the number £ being the scattering length. 
For positive scattering length there is a bound state near threshold and for 
negative scattering length a resonance ('virtual bound state'). But there is 



no finite quantum system that can give this simple model of scattering am- 
plitude: such a system can only be constructed as a limit of hamiltonians. 
Indeed this is a beautiful example of the idea of renormalization at work. 
We will construct a transfinite quantum theory of two body interactions that 
describes exactly this case of low energy scattering. 

These ideas are of interest in modern (late 1990 's) atomic physics. Tech- 
nical advances in the cooling and trapping of atoms have made it possible 
very recently to study experimentally the interaction of atoms at very low 
momenta: i.e., wavelengths large compared to the size of the atoms them- 
selves. Scattering lengths for several species have been measured. In the case 
of Rubidium and Sodium the scattering length is positive while for Lithium 
it is negative, (e.g., £ = —27.3 ± 0.8 Bohr radii for Lithium.) Interspecies 
scattering also should give a whole range of values of scattering lengths, 
both positive and negative. Indeed it is even possible to 'tune' the scatter- 
ing length to any desired value using Feshbach resonances. Thus many of 
the phenomena associated to low energy scattering have become or will soon 
become experimentally accessible. Of greatest interest is the formation of 
condensates of a large number of atoms. This would be described by a non- 
relativistic field theory (bosonic or fermionic) with a contact interaction. The 
techniques developed may eventually be useful to study phase transitions in 
such atomic condensates. So far we have a way to describe condensates in 
two dimensional systems. 

We will then pass to studying some field theoretical models. T. D. Lee 
[ ij introduced a simple model for renormalization which has been studied 
further by many people. It describes a Bose field interacting with a point- 



like source. The source itself has two possible states, which can be thought 
of as describing the internal states of a heavy particle. In the limit when 
the size of the source goes to zero, there is an ultra-violet divergence in 
this model: for example, the energy difference between the two states of the 
source is infinite. This can be removed by a renormalization procedure: the 
parameters of the model are made to depend on a cut-off ( the size of the 
source ) in such a way that the energies of physical states are finite in the 
limit as the cut-off is removed. 

An even simpler, non-relativistic, limit of the Lee model was studied in 
ref. ). The heavy particle was interpreted as the nucleon, its two states 
being the neutron and the proton. The light (but non-relativistic) boson 
was thought of as the charged pion. The energy of the 'neutron' (more 
precisely the neutron-proton mass difference) is infinite. This divergence can 
be regularized by requiring the source to have a finite size A -1 - an ultra- 
violet cutoff. Then we require a bare parameter /j,\ depend on A in such a 
way that in the limit we get a finite answer for the neutron energy. 

Now we have much better descriptions of the nucleon pion system. The 
static source model is thought of these days as a historical curiosity. But 
this kind of model has, been very valuable as a proving ground for new 
approaches to renormalization. The most spectacular example has been the 
work of Wilson, who also perfected his ideas on a variant of the static source 
model ||. A later (one dimensional) version of the static source model, the 
Kondo model, has become a classic example of renormalization. We will also 
use the static source model and (later on, the Kondo problem) to test our 
ideas on renormalization. 
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The simplest example of such a field theory we will study is a system of 
(non-relativistic) fermions interacting with a point source. This fermionic 
variant of the Lee model also has a static source with two states: one of 
these states is now a fermion and the other a boson. We will show that the 
ultra-violet divergences can be removed by a nonperturbative renormalization 
method. The energies are all finite and the theory in fact has a compeltely 
finite description, with no cutoffs or bare coupling constants in sight. The 
main lesson we learn is that renormalizable quantum field theories of this type 
( 'transfinite') also have interactions that arise from the boundary conditions 
of the wave function. Hence the dynamical information is not encoded in the 
hamiltonian, but in the the principal operator, $(£'), which we construct. 
The spectrum of the principal operator determines the energy levels: the 
Schrodinger equation is replaced by the equation Q(E)if) = 0. Although 
the problem is not exactly solvable, we can apply traditional approximation 
methods such as the variational principle once we have a finite form of the 
theory. We won't be able to get all the energy levels: but we will show 
that they are all finite and determine the ground state and first excited state 
energies. 

We will then apply our methods to the bosonic static source model (the 
original Lee model). This case is more subtle. In fact the original analysis 
of the Lee model was incomplete: the divergences were shown to be removed 
only in states which contain at most two bosons. The possibility remained 
that the hamiltonian of the system is unbounded below, when the number of 
bosons is more than two. This is not a mere technicality: there are systems 
in which an analogous renormalization model still leaves behind divergent 
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energies for multi-particle states |7J. ( We will give such an example ourselves 
later on in this paper.) We will show that the non-relativistic Lee model has 
energies bounded below in each sector with a fixed number of bosons ( this 
number is a conserved quantity of this model.) Moreover we will make a 
variational estimate of the ground state energy of the Lee model in the limit 
of a large number of bosons. Thus it will be established that the Lee model 
is indeed free of divergences. 

Next, we will apply these ideas to the case of a system of non-relativistic 
bosons in two dimensions, interacting though a contact interaction: non- 
relativistic A0 4 theory in 2 + 1 dimensions. This model has been studied in 
various guises by other authors as well ||, |, [L(J, but our approach is some- 
what different. We will obtain a closed form for the principal operator after 
renormalization. We will then show how to solve the many body problem 
in the mean field approximation. In fact we will obtain a solution for the 
wavefunction of the bosonic condensate in this aproximation, as well as its 
energy. 

2 Scattering at Low Momentum and Renor- 
malization 

In the limit of small momentum we should expect that the scattering ampli- 
tude of two atoms does not depend on the details of the atomic form-factors: 
the 'shape' of the atoms should not matter at wavelengths much larger than 
the atoms. This is a situation to which the philosophy of renormalization 
applies perfectly: the cut-off is the size of the atom: the independence of 
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the atomic scattering amplitude on the shape of the atom arises from the 
independence on the regularization scheme of the renormalized system. We 
will now show in some detail how to perform this renormalization. 

Elastic scattering |TTJ of a particle by a heavy target is described by 
the scattering amplitude f(k, n, n'). ( Of course the case of two-parcticle 
scattering can be reduced to this case by passing to the center of mass frame.) 
It is defined in terms of the asymptotic form of the wavefunction as r — > oo, 

^(k,x)~e^ n - n ' + /(fc,n,n') — • (2) 

r 

Here, k = /en is the momentum of the incoming wave and n' = - is the 
direction of the outgoing wave. Conservation of probability requires the 
"unitarity condition" on the scattering amplitude: 

/(*, n, n') - /* (fc, n\ n) = 2ik j f(k, n, n")f*(k, n\ n")^. (3) 

Let us also recall the formula for the total scattering cross-section: 

a(k)=J\f(k,n,n')\ 2 dn n ,. (4) 

If a wave is scattered by a target whose size is small compared with the 
wavelength, only the partial wave with zero angular momentum will scatter: 
only this partial wave has an appreciable probability of interacting with the 
target. Thus the scattering amplitude will become independent of angles 
in this limit, even if the target is not spherically symmetric. The unitarity 
condition will then require that 

lmf(k) = k\f(k)\ 2 . (5) 
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In other words 

lm W) = - k (6) 

Thus, in the limit of vanishing momentum the scattering amplitude ap- 
proaches a real constant (the unitarity condition requires the imaginary part 
to be O(ifc)): 

limf(fc,n,n') = -e (7) 

fe— +o 

The quantity £ is called the "scattering length". It can be either positive 
or negative. It is stated in some textbooks that a positive scattering length 
corresponds to a generally repulsive interaction and a negative one to an 
attractive interaction. But this is only true of interactions that are weak 
enough to be treated in the Born expansion. We will in fact see that a 
positive scattering length can even lead to a bound state. 

Now the imaginary part of jhr must be O(k) to satisfy the unitary con- 
dition. In the limit of small but nonzero momentum we thus have, 

wr~ ik -\ (8) 

or, 

'(*■"•»') = "FT*- 0» 

In fact the scattering amplitude of any target whose size is small compared 
to the wavelength is described asymptotically by this formula: even if the 
target has no special property such as spherical symmetry. We will call this 
the 'simple pole model' for low momentum scattering. (Due to Wigner(1933) 
and Bethe and Peierls (1935).) The scattering cross-section of the target will 
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be, in this model, 



a{k)=Alxe TV¥e- (10) 



This gives the geometrical meaning of the scattering length: at low momenta 
the target will appear to be a hard sphere of radius £. 
The scattering amplitude has a simple pole at 

*4 (11) 

This pole is in the upper half of the complex £;-plane for £ > 0; it then 
corresponds to a bound state of binding energy -^pz- If £ is negative, the 
pole does not describe a bound state. Poles of f(k) in the lower half of 
the /c-plane correspond to solutions of the Schodinger equation that grow at 
infinity: it corresponds to what is called a 'virtual level'. ( See [II] section 
133.) 

This simple picture for the scattering at low momenta however has an 
important peculiarity: there is no hamiltonian of the form 

ff =fl + V W (12) 

with a potential function V(x) which can reproduce such a scattering ampli- 
tude exactly. If there were such a potential, it would have zero range and 
infinite height (for positive £) or depth ( for negative £.) V(x) cannot be 
described even by a familiar distribution such as a delta function: the delta 
function potential in three dimensions either has vanishing scattering ampli- 
tude (repulsive case) or has no well-defined ground state (attractive case). 

Instead we have to view the hamiltonian as arising from a limit of poten- 
tials, with widths tending to zero; the heights have to be carefully adjusted 
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as a function of the widths in this limit in order to get a non- vanishing scat- 
tering amplitude or a finite ground state energy. This process of obtaining 
a finite scattering amplitude from the limit of a sequence of potentials with 
finite width is reminiscent of the renormalization program of quantum field 
theory; the width is the 'short- distance cut-off' and the height (or depth) of 
the potential the 'bare coupling constant'. The main difference from conven- 
tional quantum field theory is that the renormalization has to be carried out 
non-perturbatively, since we expect to recover a bound state in some cases. 
The answer will be independent of the details of the limiting process. 

We will now show that there is in fact a perfectly well-defined quantum 
theory with the simple pole model above as the scattering amplitude. Its 
hamiltonian is the same as that of the free particle as far its action on smooth 
wavefunctions in position space is concerned. The interactions are encoded 
into the boundary conditions of the wavefunction at short distances ( or at 
infinity in momentum space). Thus the hamiltonian is practically useless 
as a tool in studying this system. We will instead obtain a formula for 
the resolvent of the regularized hamiltonian and take its limit as the cutoff 
is removed. It will turn out that it is determined in terms of a function 
4>(E) (the 'principal function'): we get a version of the Krein formula for 
resolvents. The zeros of this principal function give the point spectrum of 
the renomalized system. 

When we have more than two particles, we can still preform the renormal- 
ization as before and get a formula for the resolvent. However, this time the 
principal function is replaced by an operator. In the case of the three body 
problem, we can study the spectrum of this operator. If the dimension of 
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space is two, we show that this problem is well-posed and has a well-defined 
ground state energy. In the case of three-dimensional three body problem, 
the ground state energy of the renormalized theory still diverges: there are 
further renormalizations necessary. Thus even if we can obtain a Krein for- 
mula and a principal operator, we still need to show that the spectrum is 
bounded below in order to have a well-defined theory. This is why we take 
pains to establish this lower bound in the case of some quantum field theories. 

2.1 Renormalized Resolvent 

Consider a pair of particles with an attractive short range interaction. After 
separating out the center of mass variable, we can reduce this to the scattering 
of a particle of mass m (equal to the reduced mass of the pair) against an 
immovable 'target' representing the interaction between the particles. We 
are interested in the limit as the inverse range of the interaction A is very 
large compared to the momentum of the particle. 

This can be modelled by the Hamiltonian operator, in momentum space, 

2 

HaHp) = f^b) - <7(A)pa(p) J PA(qMq)idq]. (13) 

Here, p\(p) is a function that is equal to one near the origin and falls off 
rapidly at infinity. For example, 

Pa(p) = 6(\p\ < A) (14) 

would be a typical choice. The function g(A) will be picked later in such a 
way that the scattering amplitude has a limit as A — ► oo. 

Here, the interaction is represented by a separable kernel which makes 
the calculations simple. If instead we choose a potential of range -^ we 
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will get similar answers in the limit A — »■ oo but the calculations are more 
complicated. 

Consider the inhomogenous equation 

2 

(|- - e)Mp) - 9Wpa(p) J PA(q)Mq)[dq] = x(p)- (15) 

Then 

MP) = -^^ + -^^PAiP) (16) 

1m 2m 

where A^ is 

A A = g(A)J P/v (p)Mp)[dp}- (17) 

We can put the expression for ijj\(q) into this equation for A A to get, 



A; 



9-\A) ~ ' PliP) 



Z^-E 



Ira ' Ira 

We now choose g(A) such that 



PA(P)#^H (18) 



^(A)-/#^H (19) 

1m 



has a limit as A — ► oo: 



g -\K) = J p\(p)^-— E [d v \ (20) 

Ira 1m 

for some real constant \x. 

This number /^( which has the dimension of momentum) is the true phys- 
ical parameter which describes the strength of the interaction: it remains 
meaningful even in the limit as A goes to infinity. If the 'bare coupling con- 
stant' g(A) is eliminated in favor of //, all the divergences will dissappear: this 
is the essence of renormalization. This kind of replacement of a (divergent) 



coupling constant by a momentum scale is quite common in renormalization 
theory. It is sometimes called 'dimensional transmutation'. 

Then A\ will have a limit as A — ► oo, and so will the solution ip(p) = 
lim A ^oo ^a(p) : 



i>(p) 



x(p) 

2m 



1 



1 



E <t>(»,E)g--E 



x(q) 

2m 



E 



[dq]. 



(21) 



Here, 



<K»,E) 



[dp] 
2m 



tJ.iL pi e 

2m 2m 2m 

" 2 + < 



[dp] 



p- 



p 2 + jj, 2 



PZ 



-2mE) 
-2mE 



f22) 



This is a convergent integral. Note that the limiting solution is independent 
of the choice of p\(p). 

Now, the resolvent kernel is given by the formula, 



ijj(p) = / R(E;p,q)x(q)[dp\. 



(23) 



Thus we get, 



mp , k) = W<<z^> + 1 



(24) 



This is called the 'Krein formula' for the resolvent. We can regard our 
limiting system as defined by this formula for the resolvent. The hamiltonian 
from which this follows is, as a differential operator, the same as that of the 
free particle. The interaction is encoded into the boundary conditions at the 
origin. The resolvent, being a Green's function, encodes the information of 
these boundary conditions as well. The function ,, l E ~, gives thus a convenient 
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description of the interaction. We will see that this is (upto a constant) the 
scattering amplitude. 

Thus the interactions of the theory are described by the function </)(fi, E). 
Apart from it, all the terms in the formula for the resolvent just involve the 
free theory. We will see that with more than two particles, we still have 
a similar formula, but the real valued function 0(/x, E) is replaced by an 
operator valued function, the 'principal operator'. This principal operator 
acts on a reduced Hilbert space which for the simple case ofthis section is 
one dimensional. That is why we just have a function 4>{[x, E) rather than 
an operator. 

2.2 The Krein Formula and Boundary Conditions 

The resolvent of a differential operator contains the information on its bound- 
ary conditions as well. It is interesting to make explicit these boundary con- 
ditions implied by the Krein formula for the hamiltonian. These boundary 
conditions determine the domain of the hamiltonian thought of as a self- 
adjoint unbounded operator in L 2 (R 3 ). 

The domain of the hamiltonian is the range of its resolvent: the set of 
functions in momentum space that can be written as 

J R(E;p,k) X (k)[dk] (25) 

for square integrable x{k)- Thus for the free particle with resolvent 

R (E;p,k)= K J ^_ E } (26) 

1m 
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we have the set 



X(P) |. .. r2 , r >^ 



W(p) = ^rr^ lx G L (i? )} - (27) 

2m 

In other words, the domain of the free hamiltonian consists of wave functions 
ip(p) for which both if) and p 2 ip are square integrable. In position space this 
means that the domain of the free hamiltonian is the set of wavefunctions 
that are square integrable and have square integrable second derivatives. 

Now let us ask how the form of the resolvent changes if we change the 
domain of the hamiltonian. The difference of the resolvent from its value for 
a free particle is ( in position space) 

R{E;x,y)-R (E;x,y) (28) 

a homogenous solution of the Schrodinger equation: two different Greens 
functions for the same differential equation differ by a homogenous solution. 
The simplest possibility is to change the behavior of the resolvent at one point 
( say the origin) by putting there a pointlike scatterer. ( More complicated 
modifications are also allowed mathematically, but are not as interesting ). 
Then we must have, for d — 3, 

iy/(2mE)\x\ -iy/(2mE)\y\ 

R(E; x, y) = R (E; x, y) + C(E) (29) 

\x\ \y\ 

for some constant C(E). We require that the additional term corespond 
to outgoing waves at spatial infinity, a physical requirement. The formula 
for arbitrary d is similar and involves Hankel functions. This means the 
wavefunction can blow up at the origin as tA=t', (or log jx| when d — 2). 
Note that such a singularity is still square integrable. 
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In momentum space this is of the form 

R(E; P , k) J 2 ^- k > + B( E) j / l . (30) 

2m 2m 2m 

The quantity B{E) is determined by the condition that this operator be in 
fact a resolvent: 

R(E;p,k)-R(E';p,k) 



E-E 1 
After some calculations we get 

B~\E) - B-\E') r 1 1 



f R(E;p,k')R(E',k',k)[dk}. (31) 



2m 2m 

Or, 

fi- 1 ^) - fi-^F) = J[dk}{^-^ - ^-^} (33) 

2m 2m 

But this is precisely what we got from renormalization: our quantity B(E) 
is just the inverse of the <p(E) we had previously. 

Thus we see that the renormalized hamiltonian is just the free hamiltonian 
with a modified boundary condition on the wavefunctions at the origin of 
position space. 

2.3 The Scattering Amplitude 

Consider the Schrdinger equation for a free particle: 

_Lv 2 ^(x) =Exj){x). (34) 

Suppose that E > 0. If we require that the wavefunction be continuos 
everywhere in space, we have the usual plane wave solution: 

Ce ik - X (35) 
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with E = Jf—. In momentum space this corresponds to a wavefunction 

(2n) d 5(p-q). (36) 

But this is not the only solution if we allow the solution to blow up at one 
point (say the origin). This would mean that there is a static scatterer 
of zero size sitting at the origin. We would still require the Schrodinger 
equation to hold away from the origin. Hence the solution would have to 
differ from the plane wave by a multiple of a homogenous solution of the 
differential equation, one that may diverge at the origin. This homogenous 
solution represents the scattering by the particle at the origin; hence we 
should require that it become an outgoing wave at infinity: 

p ikr 

4>(x)=e ik - x + f(k) — . (37) 

r 

Here the constant f(k) has the meaning of the scattering amplitude. In 
momentum space this becomes 

2ir 1 

^(p) = (2n) d S(p -k) + f(k)-- 2 - 2 (38) 

m -2 * ze 

2m 2m 

it being understood as usual that e — ► 0+. 

We can now compare this with what we get from the formula for the 
resolvent and see that the extra term in the resolvent is precisely of this 
form: 



V(p) = (2ir) d 8(p -k)+ - k2 - 2 - 2 (39) 



1 

This is in line with our argument that the Krein formula for the resolvent de- 
scribes a boundary condition on the wavefunction at the origin: it is allowed 
to blow up. 
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Thus we have a scattering amplitude that is independent of angles: 



/;/ 



Evaluating the integral, 



/(*) = n _ lt . fc2 . - ■ (4°) 



<Kn,E) = — 



^/(—2mE — ie) — fj, 



(41) 



so that the scattering amplitude is 



/(fc > = ^Tk (42) 

The sign is fixed by the rule lim^ol - 1 — ie) 1 / 2 = —i. 

Thus we find exactly the 'simple pole model' for the scattering amplitude 
for low momenta that we had in the last section, with scattering length £ = -. 

Later on we will study the many body problem of particles with such 
contact interactions. We will give a complete theory only in the simpler 
case of two dimensional space. Problems with the extension to the three 
dimensional case will be described as well. 

3 The Fermions with a Static Source 

Since our fermionic variant of the Lee model is simpler we will describe it 
first. 

Let ip(p),i>'(p) be the creation-annihilation operators for a fermion field 
in d+1-dimensional space-time. They are represented on the Fock space T 
built from the vacuum |0 >: 

m P ),^(q)} + = (2n) d 8(p-q), [tf (p) , tf (q)]+ = = [^{p),^(q)] + (43) 
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and 

^(p)|0>=0. (44) 

We are mainly interested in the case d = 3 but our method applies for any 
d<A. 

On JF <g> C 2 , define the hamiltonian 

H A = H + H 1A (45) 

where, 

#o = J[dp]i/>*(p)tl>(p)w(p) (46) 

and 

#ia = A*a^^ +gJ[dp][pA(p)i>(p)a- + /i.e.]. (47) 

There is no loss in assuming that the coupling constant g > since any phase 
in g can be absorbed into a redefinition of the field ip(p). The quantity /ia is 
a 'bare coupling constant' whose dependence on A will be determined later 
by renormalization. 

The dispersion relation of the fermion field is chosen to be non-relativistic: 

2 

uo p = m + ^-. (48) 

2m 

The Pauli matrices cr±,(T3 act on C 2 in the usual way. There is a U(l) 
symmetry (rather like isopin) which leads to the conserved quantity 

Q= l -^ + J[dpWi>{p). (49) 

Clearly Q > 0. 

This model describes the interaction of the fermions with some heavy 
particles sitting at the origin. The function Pa(p) is a sort of form- factor: it 
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describes the the internal structure of these heavy particles. For us it will 
serve the purpose of an ultra-violet regulator, eg. we may choose 

Pa(p) = 0(\p\ < A) (50) 

for some momentum cut-off A. In the limit A — ► oo we will have point-like 
heavy particles. But in this limit there is a UV divergence in the theory. 
This divergence can be removed by renormalizing the constant /j,\. The final 
answer will not depend on the choice of the form- factor p^: any function 
which is equal to one at the origin and falls of faster than any power at 
infinity will give the same answers. This is part of the universality of the 
transfinite theory: its independence on the regularization scheme. 

Our theory can be thought of as a model for the interaction of the top 
and charm quarks with the Higgs boson: the heavy boson is the Higgs boson, 
the heavy fermion the top quark and light fermion the charm quark. Our 
model is an approximation where m t , ran >> m c — "m but \m t — ra#| « va c . 
In this limit we would have the interaction of non-relativistic b-quarks with 
either a Higgs boson or a top quark; the parameter /i describes the mass 
difference betwen the two heavy particles. 

It is certainly true that that m t » m c and it is almost certain that 
itih » tti c . It would be unusual for the Higgs boson and the top quark to 
have almost the same mass, as we are asuming. Also, the Higgs-top-charm 
coupling is very small in the standard model (induced by higher loop effects) 
but is not so small in some variants of the standard model (for example with 
super symmetry). Thus our model describes a somewhat unusual possibility 
in particle physics: but this has not yet been ruled out experimentally. In 
any case we can use this system as a physical picture that guides our mathe- 
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matical analysis of this renormalization problem, much like the nucleon-pion 
system in the Lee model. More realistic cases can be studied later, once the 
principles are established in this simple case. 

In some extensions of the standard model (eg. the minimal supersym- 
metric exension) there are charged Higgs particles. In these models, a decay 
of the top quark into a Higgs boson and a bottom quark is possible. Again, 
if the Higgs and the top quark happen to be almost degenerate, our model 
will describe this system. 

The hamiltonian can be expressed as a 2 x 2 block split according to C 2 : 

H A = . (51) 

\9 J[dp\pA{PW(p) #o + Ma / 

We will now construct the resolvent of this hamiltonian using some identities 

given in the appendix. Then we will show that the limit as A — > oo of 

the resolvent exists: this is renormalization. After that we will study the 

spectrum of the renormalized theory. 

In the appendix, we work out some formulae for inverting operators split 

into 2x2 blocks as above. In the notation used there, the resolvent iof the 

regularized hamiltonian s 

«»w = jjlb=l- '-I (52) 

where, 

1 




a 



^77 F 6t$ A(£) -1 &77 F> ^ = "77 F^aO^) -1 (53) 



H $ — E Hq — E Hq — E Hg — E 

and 

6 = $ A {E)- 1 , b = gJ[dp]p A (p)i;(p). (54) 
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We express everything in terms of § A (E) since we will see that all the di- 
vergences are removed once we have a proper definition for Q A (E). Once we 
know $, the resolvent is given by the above explicit formula. 
Indeed, we have, 

$ A (£) = H - E + p A - g 2 J[dp][dq]p A (p)p A (q)^(p)j^-^^(q). (55) 

Now we normal-order the last term; i.e., we reorder the operators so that the 
creation operators stand to the left of the annihilation operators. This can 
be done using the identities 

H ^(q) = tf(q)[H +u(q)], Hp)^(q) = -tf(<l)1>(p)+(2*) d 8(p-<l), (56) 
Thus 

$ A (E) = Ho-E + g 2 f[dp][dq]p A (p)p A (q)^(p) * *ft(g) 

J H + uj{p) + uj{q) - E 

+»*-!> 2 JWpI(p) Ho + w \ p) _ e (57) 



*0 

Suppose we now choose 



Pa{p) 2 



, 2 //W [dp ft . a \ 

Here yU is a parameter independent of A. Moreover we will require that m > p 
for simplicity. (This is not essential: it guarantees that both the heavy states 
are stable against decay.) 

The point of this choice of p A is that it cancels the divergent part of 
$ A (£): 

* A (E) = H -E + g 2 ([dpdqWip) - - pA ^ff Mq) 

J H + uj{p) + u>{q) - E 
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-9 2 / [dp}pA(py[ H - - — ] + n (59) 

J H + uj{p) - E u(p) - fi 

The integrand of the last term behaves for large \p\ like / I P L '■ it is 

2 

convergent if d < 4, since u = m+ %—. Thus we can take the limit as A — *■ oo 
keeping \x fixed at a positive value less than m. Then, the limiting operator 

$(£) = H -E + g 2 f[dpdq}^(p) \ -<%) 

-9 2 f[dp][ „ , ) , r- - — ^ ] +/i (60) 

exists. This is the renormalized form of the principal operator. W can just 
put this into the earlier formula to get a formula for the resolvent. 



3.1 The Principal Operator 

To get more explicit expressions we specialize to the case d — 3. (In fact the 
following arguments apply to any d < 4.) Evaluating the integral, 

$(E) = H -E + ^ + 47T 2 g 2 (2m)2y(H + m-E)-^(m-^)] 



+g [dpdq}^(p)—— -— — -i/j(q) (61) 

J H + uj{p) + uj{q) - E 

This principal operator can be written as the sum of a 'kinetic' (single par- 
ticle) term 

K(E) = H -E + n + ATi 2 g 2 (2m)iy(H Q +m-E)- y/(m - »)] (62) 

and an 'interaction' term 

U(E) = g 2 f[dpdq]^(p) - ■ . .\ . . Mq). (63) 



H + uj(p) + w(q) -E 
The fermion number operator N = f[dp]ifi< (p)?p(p) commutes with $(£'). 
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The principal operator determines the dynamics of the theory. The re- 
solvent of the hamiltonian of the renormalized theory is 

(a{E) I3\E) 
R{E) = | (64) 

U{E) $(£)-! ' 



where 



° (E) = H^E + W L E h^W) JW^H^E < 65) 



and 



(3(E) = --i-r f[dp]cf>(p)- ' 



$(E)J l ^^'H - E' 
This can be thought of as a sort of Krein formula for the resolvent of our 
field theory. 

Indirectly it defines the hamiltonian H as the operator for which R(E) = 
H l _ E is the resolvent. But if we were to think of this renormalized hamiltonian 
directly as an operator on the fermionic Fock space, it would appear to be 
just the free fermion operator. For example, consider the space of states 

1/ >= Jfipi ■ ■ -PnWiPl) ■ ■■^(Pn)[d P l ■ ■■d Pn ]\0 > (67) 

with smooth f(p±, ■ ■ -p n ) decreasing at infinity faster than any polynomial. 
They correspond to fermionic wavefunctions that are smooth in position 
space that also fall of faster than any polynomial. On these 'nice' states, 
the hamiltonian H of the transfinite theory is the same the free hamiltonian: 

H\f >= /[w(pi) + • • -u;(p n )]/(pi • • -PnWipi) ■ ■■^(Pn)\0 > [dpi ■ ■ ■ d Pn ]. 

(68) 

But this formula does not uniquely define if as a self-adjoint operator. It 

has many self-adjoint extensions, corresponding to different sets of boundary 
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conditions as the momenta go to infinity. One of them is the free hamiltonian 
Ho, but there are other extensions too which define interacting theories. The 
interactions are hidden in the domain of definition of this unbounded opera- 
tor: in other words in the boundary conditions of the fermion wavefunctions 
in the limit of large momentum. It is quite awkward to think this way: for 
example it would be difficult to calculate the spectrum of the renormalized 
theory. 

But the same information is contained in a more explicit form in the above 
Krein formula for the resolvent. The resolvent is a sort of Green's function 
which therefore contains also the information on the boundary conditions. 
Another way to understand it is that (away from the spectrum) the resolvent 
is a bounded operator, so there is no need to specify its domain. 

Transfinite field theories are thus theories in which the interaction takes 
place at infinity in momentum space. We are still able to give a sensible 
description of these theories. 

3.2 The Spectrum of the Transfinite Theory 

We saw that there is an explicit formula for the resolvent of the hamiltonian 
in terms of the inverse of the Principal operator. Thus, all the dynamical 
information about the theory is contained in its principal operator. For 
example the discrete spectrum of the hamiltonian correspond to the poles of 
the resolvent. There are no poles in H l _ E : its spectrum is purely continuous. 
Thus the poles must arise from those of $(i?) _1 ; i.e., roots of the equation 

$(E)\u>=0 (69) 
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This equation now plays the role of the Schrodinger eigenvalue problem in the 
transfmite theory. The residue of the pole of the resolvent is the projection 
operator to the corresponding eignspace of H. Thus the eigenvector of H 
corresponding to a root E of the principal operator is given by 

J[dp]^(p)\u>\ 

1 • (70) 

This is because the residue of a pole in the resolvent is the projection operator 
to the eigenspace with that eigenvalue. We can read off this residue and see 
that it is the projection to the above state, once a root of the equation 
$(E)\u >= is found. 

An example is the vacuum state: 

^(E)\0>=(-E + fi + g 2 (2m)2y(m-E)-^/(m~i2)})\0> . (71) 

There is a root when 

E = n. (72) 

Since // < m, it is not in the spectrum of Hq. The corresponding eigenvector 
of the hamiltonian is 

a-f[d P W(p)\0> 



Ho-fi 



(73) 
> 

It contains a fermion in the first entry, so it is not the vacuum in the 
whole Hilbert space T ®C 2 . In spite of this, it is the state of lowest energy- 
the ground state- when /x < 0. (Proof is in the next section). When m > 
2 This root is assumed not be in the spectrum of Hq. In this case it would be be 
embedded in the continuous spectrum of H and should be unstable. We will deal with the 
continuous spectrum of H later. 
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fi > the above state is the first excited state. ( The ground state is just the 

/|0>\ 
'vacuum' state .) When fj, > m, we don't get an eigenstate of the 

V o J 

hamiltonian since this state is then unstable. (The ground state in this case 

/|0> 

is also 

v ° 

More generally, the spectrum of the hamiltonian is the set of values of 
E at which the resolvent either does not exist (discrete spectrum) or exists 
but is unbounded (continuous spectrum). Thus the continuous spectrum will 
be that of H plus the values of E at which $(£?) does not have a bounded 
inverse. 

3.3 Proof That The Ground State Energy is Finite 

In order to see that we have exorcised all the infinities, we must show not 
only that the resolvent of the hamiltonian exists in the limit as A — ► oo but 
also that the ground state energy is finite. This is nontrivial to prove since 
we will see that there are many theories where even after a renormalization 
there are further divergences which make the spectrum not bounded below. 

We will estimate the norm ||<I>(.E) _1 || for the case E < /x; our aim is to 
show that this is finite. Then the ground state energy is either zero (when 
H > 0) or is equal to /i. 

It is sufficient to consider the sector with the fermion number held fixed 
at some value n, since <$>(E) preserves this number. 

Recall that $(£) = K(E) + U(E) with 

K(E) = H -E + /j + ATi 2 g 2 {2m)iy{H Q + m - E) - y/(m - (i)] (74) 
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and 



U(E) = i I [dpdqWip)li — { -l—^-^ m . (75) 



Moreover, 

K{E)>nm + {n-E), U(E)>0. (76) 

The inequalities become equalities when the fermion number is zero. 

Now write, 

$(E) = K(E)*[1 + U(E)]K(E)*, U(E) = K(E)-*U(E)K(E)~*. (77) 

Since U(E) > 0, 

ME)- l \\ < \\k(e)-'\\ B + um-'w 

< [nm + (// — -E)]~ 

< -^-p- (78) 
IX — E 

This shows that $(£') has spectrum bounded below by /i. 

This simple method is not sufficient to determine the ground state in the 
bosonic case. We will have to supplement it with a sort of mean field theory. 

4 The Lee Model 

The Lee model has a charged bosonic field satisfying: 

[0(p), 0%)] = (^YHp ~ Q), W{P), <PKq)) = = [0(p), 0(?)]- (79) 
The bosonic Fock space £> is built from the vacuum in the usual way: 

0(p)|O >= 0. (80) 
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Again, the complete Hilbert space of the system is B ® C 2 . 
The hamiltonian is, again, 

H A = H + H 1A (81) 

with 

H = J[dp\4>\p)4>{p)uj p (82) 

and 

#ia = Ha — g - ^ + 9 J [dp][pA(p)<f>(p)<r- + h.c.]. (83) 

Some of the methods are the same as in the fermionic case; we will then omit 
the details. But the bosonic case requires some new analysis as well. 

The divergences are removed as before by normal ordering and choosing 
the bare mass difference to be 

2 f pA(p) 2 [dp] 
J w(p) - n 

In the limit the resolvent has the form: 

1 fa(E) (3\E) \ 

R(E) = lim tt - =\ ; (85) 

«<*> = H^E + ^I^^mI MM B^E' t m 

m = ~S^^^ (87) 

and finally 

$(£) = H -E + /j + 47r 2 g 2 (2m)ly(H + m-E)-^/(m-^} 

-g 2 [dpdq]${p)—— — -— — — -0(g). 

J H + uj{p) + Lu(q) - E 



35 



Other than the sign of the last term the answer is essentially the same as 

before. But this sign makes an important difference. 

/|0>\ 
The proton is the state and it is an eigenstate of the renormal- 

v o ; 

ized hamiltonian with eigenvalue zero. The neutron is an eigenstate of the 

hamiltonian, 

^L_/[dp]0t(p)|O>\ (gg) 

It corresponds the root of 

$(E)\0 >= (-E + ii + g 2 (2m)%[^(m-E)- y/(m-fi)])\0 > (90) 

with E — \x. Thus the parameter /i is the renomalized value of the neutron- 
proton mass difference. So it is reasonable to assume m > \i > 0. 

But is the proton the state of least energy? Could there be states which 
contain many bosons that have a lower energy?. Is there even a ground state? 
These questions can be answered by studying the principal operator $(£"). 

4.1 Lower Bound for the Ground State Energy 

We will estimate the norm ||$(£) -1 || for the case E < \i. It is sufficient to 
consider the sector with the number of bosons held fixed at n, since <&(E) 
preserves this number. 

Let us define the "kinetic" part of $(£") to be: 

K(E) =H -E + ii + 4n 2 g 2 (2m)%W(Ho + m-E)-^/(m-n)]. (91) 

Clearly, 

K(E)>nm + {n-E). (92) 
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Define -U(E) to be the "potential" part of $(E): 



$(E) = K(E)-U{E). 



(93) 



Notice that <&(E) is now the difference of two positive oeprators rather than 
the sum as in the fermionic case. This is what makes the bosonic case more 
complicated. 
As before, 

$(E) = K(E)*[1 - U{E)]K(E)^ (94) 



where 



Also, 



i.e., 



U(E) = K(E)~2U(E)K(E)~s. 



I^r 1 !!^!!^)- 1 !! \\[i-u(E)} 



-ii 



imEy'WKlnm+ifi-E^Wil-UiE)}- 1 ] 
Now, by explicit calculation, 

U(E)=g 2 [[dpdq]<j>t(p) [ 



K(E - w p )3(flo + uj p + u q - E)K(E - u^ 



(95) 



(96) 



(97) 



(98) 



Thus, (remembering that inside the square bracket the boson number is 

n-1), 



U{E) < g 2 [dpdq](j)\p)[(n - l)m + \x- E + u p ]~* 
[(n - l)m + ll> p + uj q - E]~ 
[(n - \)m + fi- E + u q ]~*<f>(q). 



(99) 
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Also, in the sector with n bosons, 

(/> (p)u(p,q)<j)(q)[dpdq]\\ < n 



\u(p,q)\ 2 [dpdq] 



(100) 



Combining these, we get 



\\U(E)\\ 2 < g A n 2 J [dpdq][{n - l)m + [1 - E + Up}' 1 
[(n - l)m + uj p + u q - E]~ 
[(n - l)ra + n - E + cu q ]~ . 



(101) 



Now we put in u(p) = ^— + m.In the middle factor we can replace an m by 
a \x since m > ji: 

2 2 2 2 

[(n-l)m + 2m+^ + ^-- E}~ 2 < [nm + // + ^- + ^- - E}~ 2 . (102) 
2m 2m 2m 2m 

This makes all the constants in the denominators the same so that we can 
scale them out to get 

(2m) 3 



4„2 



\U(E)\r < 9*n 



[nm + fi — E] 
[dpdq] 



;i+p 2 )(l + g 2 )(l+p 2 + ? 2 ) 2 
The integral is convergent. Working our way back to the beginning, 

1 1 



(103) 



mE)\r< 



where 



y/{nm + a — E) 



C 2 



\J{nm + fi — E) — g 2 n(2m) 2 C 



[dpdq] 



(104) 



:i+p 2 )(l + q 2 )(l+p 2 + q 2 y (1 ° 5) 

If E is to be an eigenvalue of energy, it must be big enough to make the 
denominators on the r.h.s. vanish; otherwise, $(£') _1 would remain bounded. 
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This means there is a lower bound on all eigenvalues, which really is the same 
as a lower bound on the ground state energy: 

E gr > (nm + fi)- n 2 g 4 (2m) 3 C 2 . (106) 

Thus we see that in each sector with a fixed number of bosons, there is a 
ground state. However, there is still the possibility that the ground state 
energy diverges as n grows to infinity. In this limit we should expect all the 
bosons to settle into the same state. The ground state of the free theory 
has all the bosons in the zero momentum state. In the interacting theory, 
the ground state could in general be something quite different. For example, 
the bosons might settle into a state which is concentrated at the origin. 
The difference between that bosonic number density and its free field value 
is called the 'boson condensate'. Whether such a non-zero condensate of 
bosons forms cannot be settled by the present analysis: we need to study the 
limit as n — ► oo. 

4.2 The Large n Limit of the Lee Model 

In the limit that the number of bosons becomes large, we should be able to 
use mean field theory. We would expect all the bosons to occupy the same 
state u(p). This state is normalized so that the occupation number is n: 

\\u\\ 2 = f\u(p)\ 2 [dp] =n. (107) 



In the limit n — *■ oo, operators can be approximated by their expectation 
values in this state: a kind of mean field theory. They will then become 
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functions on the space of such states, the complex projective space 3 of L 2 (R 3 ). 
Thus our principal operator becomes the principal function 

$(E,u) = h (u)-E + n + 4Tr 2 g 2 (2m)*W(h (u)+m-E)-^(m-n)] 

J h (u) + uj[p) + uj{q) - E 

Here, 

ho(u) = J u(p)\u(p)\ 2 [dp]. (109) 

We must solve <&(E, u) = to get E as a function of E. Then we must find 
the u that gives the smallest such E, subject to the constraint on the norm 

INI 2 = n - 

It is convenient to reexpress the problem in terms of some new variables. 
The principal function depends on E only through the combination 



A = / u(p)\u(p)\ 2 [dp] - E. (110) 

So define, 

f(X,u) = X + i2+(27rg) 2 (2m)l[^(X + m)-^(m-ii)} 

-a 2 [\dpda] M *frMg) (Hi) 

Then A is determined by the equation /(A, u) — and E by 

E = nm+ [ P 2 \u(p)\ 2 [dp] - A. (112) 

2m J 

Putting in the explicit form of to(p), 

f(X,u) = A + /i+(27r( 7 ) 2 (2m)l[ v /(A + m)- v /(m-^)] 



3 There is no physical effect if u is replaced by e u. The space of vectors of fixed length, 
modulo a phase, is complex projective space. 
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-2mg 2 [[dpdq] 2 fbMQ) (n3) 

y J L y > 2 + g 2 + 2m(A + 2m) V ; 

We now make the change of variables 

u (p) = ^/ n [2m(2m + \)]*v(y/[2m(2m + A)]p). (114) 

The change of scale of momentum will simplify the denominator in the inte- 
gral for f(X,u). The overall factor of ^Jn turns the normalization condition 
into 

r \v(p)\ 2 [dp] = 1 (115) 



so that we can separate out the n dependence. 
Now we have 

E = nm - A + n(X + 2m) f p 2 \v(p)\ 2 [dp], (116) 

and 

f(X,u) = \ + ^ + {2^gf{2m)iy{\ + m)-^{m- l x)] 



-ng 2 {2m)^[\ + 2m] / [dpdq] 2 _T 2 , V (117) 



v*(p)v(q) 
'p 2 + q 2 + l 



Imposing /(A, u) — will give an equation for A: 



[A + 2m] — {A + ^ + (2^) 2 (2m)5[ v /(A + m) - ^{m - fj,)]} 

v*(p)y(q) 

'p 2 + q 2 + l 



ng 2 (2m)l j [dpdq] \ ^ ^ . (118) 



Then A determined as a function of v. The l.h.s. can be seen to be a 
monotonically increasing function of A by writing it as 

m 
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(Recall that p, < m). Thus we can solve for A in terms of the r.h.s. as a mono- 
tonically increasing function as well. Its minimum values is —p,, attanined 
when the r.h.s. is zero. Thus, we have 

\ = -p + f 1 (nU) (120) 

where 

a = /(2 ro )l/[«4^M£l (121) 

J p z + q z + 1 

and /i is monotonically increasing. (We won't need an explicit form for /j.) 

Now put this it into the expression for energy: 

E = nm + fi + (2m - p)K + f^nU^nK - 1] (122) 

where 

K= f[dp]p 2 \v{p)\ 2 . (123) 

Now suppose we replace v(p) by v a (p) = a~^v(a~ l p). Let K(v a ) = 
K a (v), E{y a ) = E a (v) etc. Now, K a = a 2 K and U a = a 4 U: as a function of a 
both K a and C a are increasing. Then E a is also a monotonically increasing 
function of a; the minimum for £" a will occur at a = and K,U = 0. Thus 
we find that the groundstate energy is just 

E = nm + fi. (124) 

All the bosons are in the zero momentum state in the large n limit: there is no 
condensate. In other words the ground state of the bosonic field is essentially 
the same as in the free theory: the interaction is not strong enough to modify 
the ground state substantially. We will see examples later where the ground 
state is affected by the interaction. 
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5 The \($ +1NR Model 



So far we considered theories where the coupling constants did not need 
to be renormalized. In fact the UV divergence was removed by a normal 
ordering of the principal operator. Now we will study a system where the 
coupling constant needs to be renormalized: in fact one that is asymptotically 
freef^, pfl [|T3"|1 fll4]| We will discover a trick of introducing fictitious degrees of 
freedom ('angels') which will help us reduce the renormalization again to a 
normal ordering of the principal operator. 

The theory of interest is the nonrelativistic scalar field theory in two 
space dimensions, with a A0 4 interaction. This is the many body problem 
of non-relativistic bosons interacting through a 'delta function' potential. 
It is well-known that, the two body problem has an ultra-violet divergence 
which can be removed by a coupling constant renormalization. In Ref. |T5 



the 'few body problem' was studied by quantum mechanical renormalization 
methods. In this section we will renormalize this model by viewing it as a 
non-relativistic quantum field theory: in other words we will study the many 
body problem. We will carry out the analysis in part for arbitrary d; this 
will show why our renormalization method is not sufficient for d = 3. 

Define, on the Bosonic Fock space £>, the regularized hamiltonian operator 

H A = H + H 1A (125) 

where 

2 

H = J^{p)(Hp)[dp] (126) 

and 

#ia = -0(A) / [dp 1 dp 2 dp' 1 dp 2 ]p A (pi - p 2 )p\{p' 1 - p 2 ) 
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(2tt)^(pi +p 2 -A - P ' 2 )<l?(piW(p2)<t>(p'x)<l>(p'2)- (127) 

Here 

Pxip) = 0(\p\ < A) (128) 

as before. The dimensionless constant g(A) is positive which corresponds to 
an attractive interaction between the bosons. We have shown elsewhere that 



| T5| there is an ultra-violet divergence as A — > oo, which can be removed by 
renormalizing the coupling constant. Our aim is to get a manifestly finite 
expression for the resolvent of the renormalized quantum field theory. 

Now we will introduce a little trick that simplifies our problem: we will 
introduce new particles called angels that describe a bound state of a pair 
of bosons. However, they are created by operators with unusual defining 
relations: these are necessary to avoid over-counting the number of degrees 
of freedom. 

5.1 Angels 

Define operators satisfying 

X(P)XHQ) = (2*) d 6(p - q), X (p)x(q) = = X ] {p)x\q)- (129) 

Note that it is the product and not the commutator that appears here. These 
operators can be represented on the Hilbert space C © L 2 (R d ). We can 
regard xKp) as creating an entity (we will give it the somewhat whimsical 
name 'angel') out of the empty state represented by C. There can be at most 
one angel in any state: that is the meaning of the product of the creation 
operators being zero. (This is an extreme example of the exclusion statistics 
considered in some other contexts.) 
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Now consider an augmentation of the Bosonic Hilbert space, B = B(BB 
L 2 (R d ). On it define the hamiltonian 

H A = H Il + I [dp 1 dp 2 dp 3 ]p A {p 1 - p 2 ) 

t (p 1 )0 t (p 2 )x(p3)(27r)^(pi + p 2 - ps) + h.c 



+ ^n,(i30) 



Here il is the projection operator to the subspace containing no angel : 

n = J[dp] x (p)x\p), (131) 

and Hi the projection operator to the subspace with exactly one angel: 

Hi = f[dp} X \p)x(p)- (132) 



The point of introducing angels is this: the projection of the resolvent 
of the H A , to the B (the sector with no angels) is just the resolvent of the 
original bosonic system. But we will get another formula for this resolvent, 
using angels, which has a finite limit as A — »■ oo. Indeed we will see that 
the infinity is avoided by keeping the energy of the bound state of a pair of 
bosons fixed in this limit; an angel is essentially such a bound state. 

Let us split the Hilbert space according to the angel number, with a 
corresponding splitting of the operator: 

fa &t\ 

H k - MI = (133) 

\b d ) 

with 

a:B^B, tf :B®L 2 (R d )^B,d:B®L 2 (R d )^B®L 2 (R d ). (134) 

Define an operator R\(E) split in the same way: 

1 lot /?t' 

^ E) =H^wr(e J- (135) 
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We claim that 

a= H^E ] (136) 

i.e., R\ projected to B is just the resolvent of H^. To see this, we use the 

formula we obtained earlier: 

a=[a-tfd- 1 b]- 1 . (137) 

For us now [], 

a = H -E, d=——, 
tf = l{dpidp 2 ]pK{pi-P2)4 > \pi)4 > \p2)x{Pi + 'P2) (138) 



This gives 

b ] d~ l b = g(A) jldp 1 dp2dp 3 \p A (pi-p 2 )<p j {pi)<f)' l {p2)x{Pz) 
(2n) d 5( Pl +p 2 -p 3 ) 

fWMdp'slpAiti -p'2)xHPs)<l>(Pi)M)(^)Wx+P2-Pz) 

(139) 

If we use 

X(p)x\p') = (^) d s(p-P) (140) 

we will get the required result. 

But we have another formula for this resolvent: 

a = a- 1 + crVfd - barty]- 1 ^ 1 . (141) 

This will give, 



H K -E 



a 



+ -o _1 6 t $A(^)" 1 &a- 1 (142) 



We hope there will be no confusion in using the same symbol d both for an operator 
and the dimension. 
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where 



$ a(£) = ~TI\~ [ d Pi d P2 d P'i d P2]pA(Pi-P2)PA(p'i-P2) 
A J 



X (Pi +P2) 



4>(pi)<f>(p2) H _ e 4> ] {p'i)4>\p'i 



X(pl+p' 2 )(143) 



We can regard b as the operator that converts a pair of bosons into an angel. 
Then $a(-E) is a kind of effective hamiltonian in the sector with one angel 
and two fewer bosons: its zeros are energy levels of the manybody problem. 

Note that in this way of writing the resolvent of H^, the coupling constant 
appears additively! Its renormalization can be done by separating out a 
divergent constant from Q\(E). We will do this by normal ordering the 
operators in $\(E). 

Using the canonical commutation relations and 

H Q <p\p) = <p\p)H» + 1P - (144) 

we can rewrite the quantity in the square brackets above equation as 

4>\pi)4>\p2)-fi—, — rr^ — rr^ — r~v^ — r~\ — wtfaWfa) 

H + w(pi) + u>{p 2 ) + w(pi) + uj{p 2 ) - E 

+ (2n) d 5( Pl -p' 1 ) ( p\p' 2 ) ] , 6{p 2 ) 

H + w(pi) + u>{p 2 ) + uj{p 2 ) - E 

+ (2n) d 5(p 1 -p' 2 ) ( p\p' 1 ) 1 —— <P(p 2 ) 

H + w(pi) + u>(p' 2 ) + uj{p 2 ) - E 

+ (2n) d 5(p 2 -p' 1 )^(p' 2 ) ] 6( Pl ) 

H + w(pi) + cu(p 2 ) + u){pi) - E 

+ (2n) d 5(p 2 -p' 2 )^(p' 1 ) ] 6( Pl ) 

(27r) d 6(p 1 -p' 1 )(2n) d 5(p 2 -p' 2 ) ' 



H + cuip'i) + lu(p' 2 ) - E 



(2n) d 8( Pl - p' 2 )(2ir) d S(p 2 - P [) ' 



if + w(pi) + uj(p' 2 ) - E 



(145) 
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This gives, 



$a(E) = g 1 (A) - j [dp 1 dp 2 dp' 1 dp' 2 ]p A (pi - P2)p\{p' 1 - p' 2 ) 

(ft {p[ ) </> f (p' 2 ) — - T — - T — -—^ -— 

1 



+4(27r) d 5(p 1 -p' 1 )^(p 



tiV 



+2(2n) d 5(p 1 -p' 1 )(2n) d 5(p 2 -p' 2 ) 

x(p'i+p 2 )- 



H + w(pi) + w(p 2 ) + w(pa) - E 
1 



0(Pl)0(P2 
0(P2) 



flo + w(pi)+w(^)-E 



(146) 



Now we can take the limit as A — ► oo. The only divergent term is the 
last one. If we choose for g(A) the expression from the two-body problem, 



^W = ^j AiP)-^ 



(147) 



P A + /i 2 

this divergence will cancel yielding a finite expression for <fr\(E) as A — ■> oo. 
In fact, we have, 

1 * ] e(y,^o+^-^). (148) 



lim /[dg]p A (g)[dg] 



f. _i_ ii. 

2 ~ 2 



21 _i_ i£l iJ n _|_ E1±21 



E 



Here, 



?V 2 ' 2 y 



£/w 



[dp] 



4 
2^ 



p 2 + fl 2 p 2 + is 
v 2 



p 2 (p 2 + z/ 2 ) 



2 

[dp] 



V 



p 2 (p 2 + /i 2 



(149) 



We will get 

2 2 

$(£) = jldp\ x \p)Z{^,H + P --E) X {p) 
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- / [dp 1 dp 2 dp' 1 dp' 2 ]x f (pi + P2) 

H + uj{p\) + u(p' 2 ) + w(pi) + w(p2) - E 
+A{27,) d 5{p 1 -p ! 1 ) ( t ) \p , 2 ) ) - <f>(p 2 ) 

Xip'i+ti)- (150) 

Notice that the dependence on the cut-off of the previous operator is traded 
for a dependence on the renormalization scale /j,. 
The resolvent operator has the finite form 

R(E) = + 6t l h \ (i5i) 

y ' H -E H -E $(fi,E) Hq-E K ' 

which is the analogue of the Krein formula in the case of the manybody 
problem. 



5.2 The Principal Operator 

Now let us understand the spectrum of the theory in the special case d = 2. 
We will be mostly interested in eignestates with E < 0: the bound states of 
the system. In each sector with a fixed number of particles we expect the 
ground state to be of this form. ( We will make some remarks later on the 
case d — 3 where we will see that the ideas in this section will not work in 
that case.) 
Then, 



ac-r-) = ^- 2 . (152) 



,fi 2 u\ 1 , v 2 

— , — ) = — In — 

2 ' 2 ' 4tt \i 

We can rescale all the momenta by ^/(|£^|) to get 



$0u,£) 



^ln^+ W 
2n n l 
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(153) 



where 



W = — J [dp]x ] f O)log 



H + u(p) + l 



xip) 



[dpidp 2 dp[dp' 2 ]x ! (pi +P2 



1 

H + w(^) + u{p' 2 )+uo(p 2 ) + l P[ ' >2 ' 
Xip'i+A)- (154) 



+4(2n) d 5(p 1 -p' 1 )^(p' 2 ) T _ j ■ , - ■ ; - ■ , - ■ > (p : 



Thus finding the bound state energy of our manybody system amounts 
to finding the eigenvalues of W: 

W\ip >= wip, E = -^e~ 2nw (155) 

The special case of the three body problem was studied in previous papers 



fL5| . We showed not only that the ground state energy is finite, but that it 
can be estimated by a simple variational ansatz. Since we have described it 
elsewhere, we won't elaborate on this point here. 

5.3 Bosonic Condensation in Two Dimensions 

We will now discuss the case of a large number of particles: how to do 
transfinite quantum many body theory. We will only discuss the mean field 
approximation, which should be good in the limit of a large number of par- 
ticles. In the ground state, we should expect all the bosons to condense to a 
common state u(p). Thus the problem is to determine the pair of functions 
u(p), ip(p), where ip(p) being the wavefunction of the angel. The wavefunction 
u(p) describes a condensate of bosons, which breaks translation invariance. 



50 



Of course there is no Bose-Einstein condensation for free bosons in two di- 
mensions. What we will show is that with an attractive interaction of zero 
range there is in fact such a condensation. 

The expectation value of the operator W in the state \u >,ifj > becomes, 
(for large n) the 'principal function' U, 



U 



2^ / [dp\\ij(p)\ 2 \og[nh (u) + u(p) + 1] 



- / [dp 1 dp 2 dp' 1 dp' 2 \ip(p l + P2)*i/>(p'i + P 2 ) 

u* (p[)u* (p' 2 )u(pi)u(p 2 )n(n - 1) 

nh (u) + u(p' 1 ) + u(p' 2 ) + u(pi) + u{p 2 
+4 (2TT) 2 5(p 1 -p' 1 )nu*(p' 2 )u(p 2 



nh (u) + a; (pi) + u)(p 2 ) + u{p 2 



1 



where 



(156) 



(157) 



ho(u) = J \u(p)\ 2 uj(p)[dp]. 
This U is to be minimized subject to the normalization conditions 

u(p)\ 2 [dp]=Jm P )\ 2 [dp] = l. (158) 

To get more explicit answers, let us ignore all except the leading terms as 

n — ► oo: 



U = -n 



Ji>*(pi + P2)u(pi)u(p 2 )[dpidp 2 



ho(u) 



+ — logn + 0(n°). (159) 

Z7T 



Thus the ground state energy E n of such a system of n particles is given by 

2-Trn 

£^ = V— [Ci + 0(-)] (160) 

n n 



where 



^ = inf 



I \^{p)\ 2 [dp] J \u{p)\ 2 [dp] j u(p)\u{p)\ 2 [dp] 



fip*(pi + P2)u(pi)u(p 2 )[dpidp 2 



(161) 
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We will determine £ by solving the variational problem. Determination of 
the C\ takes more work which we will not carry out here. : it is not needed 
to determine the leading large n behavior 

It is more convenient to use the corresponding expressions in position 
space, with u(x) = J u(p)e ip ' x [dp] etc. 

/ \i(j(x)\ 2 dx J \u(x)\ 2 dx J \\Vu(x)\ 2 dx 



£ = inf 



Eliminating ip gives 



where, 



I[u] 



J tj)*(x)u 2 (x)dx 
£ = inf I[u] 

u 

J \u(x)\ 2 dx J \\Vu(x)\ 2 dx 



(162) 



(163) 



(164) 



/ |w(a;)| 4 da; 
We can see that this amounts to solving a nonlinear differential equation, 

which can be derived by minimizing log I [it]: 



V l u - flu + g\u\ l u = 0. 



Here, 



P 



J\Vu\ 2 d 2 x 



J \Vu\ 2 d 2 x 



(165) 



(166) 



J \u\ 2 d 2 x ' " j \u\ 4 d 2 x 
The above partial differential equation (but without the constraints on (5 
and g) has been studied in |i~6l , Theorem 6.7.25. ( We just need the special 
case of dimension two.) There is a normalizable solution for every positive (3 
and g. In fact it is enough to find the solution for a particular pair of values 
of these constants. For example, let U\ be a solution of 



V 2 ui — u\ + Iwi| 2 x2i = 0. 



(167) 
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Then, u(x) = au\(bx) will solve the general equation with 

/3 = b~ 2 , g = a 2 b-\ (168) 

The quantity I[u] is invariant under such scale transformations: I[u\\ = I[u]. 
So it is enough to solve the special case. 

Being the ground state of a many body problem, we should expect the 
solution to be real and have circular symmetry around a point . Thus we 
only have to solve an ODE, a kind of nonlinear Bessel's equation 

v"(r) + — -v(r) + —^- = v(r) (169) 

with u(x) = yjrv(r) and r = \x\. That a square integrable solution exists is 

proven^ in ref . JT^| . 

We can get a value for £ ~ 12 by numerical solution of the ODE. We also 

get this way the shape of the 'soliton' (or 'condensate' ); i.e., the wavefunction 

of the bound state of a large number of bosons with pointlike interactions. 

We plot below the solution u(r); as shown above, the scales of the u and r 

axes are not physically relevant. The wavefunction is peaked at the origin and 

decays exponentially at infinity. The ground state thus breaks translation 

invariance but not rotation invariance.f] 

5 There is a typo in equation f of ref. (l6|, p. 384; the ^v(r) term is missing. It does 

not appear to affect the rest of the argument. 

6 I thank Govind Krishnaswami for solving the ODE numerically and producing the 

graph. 
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Ground State Wave Function 




2 4 6 8 

One important conclusion of our analysis is that the magnitude of the 
ground state energy grows exponentially with the number of particles: 



E n ~ -^ee n . 



(170) 



If we had a non-singular pair-wise interaction, we would expect the magni- 
tude of the ground state energy to grow like the number of pairs; i.e., like 



n 2 . 



5.4 Vortices in Superconductors 

A limiting case of the Landau-Ginzberg theory of superconductivity gives a 
realization of the example we have been discussing. The lagrangian of this 
field theory is 

1 



\V(p\ 2 d 3 x+- J |curL4| 2 d 3 x+ / A(|0| 



2 -a 2 ) 2 . 



(171) 
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Here, </> is a complex valued function on i? 3 , A a covariant vector field de- 
scribing the gauge potential, and 

V0 = d(p - ieA(j> (f 72) 

the covariant derivative. If A > f (Type II superconductor) this Lagrangian 
has a static solution that is cylindrically symmetric around a point in a 
plane R 2 and translationally invariant in the orthogonal direction, carrying 
a unit of magnetic flux. Alontthe axis of symmetry, the field <f> vanishes. 
Two such vortices have a repulsive interaction: the energy decreases as the 
position of the zeros of <p move further apart in the plane. ( For given flux 
per unit area, the solution that minimizes energy density is a triangular 
lattice-the Abrikosov lattice.) As A — ► 1 + the force vanishes and the energy 
is independent of the position ofthe zeros. (As A — > 1~, the force vanishes 
at finite distances, but leads to a residual attractive contact interaction.) In 
fact for any integer n there is a solution carrying n units of magnetic flux, 
with zeros of (ft at any prescribed set of n points on the plane [T9| . 



Thus this vortices behave like identical free particles, located at the po- 
sition of the zeros of the field (p. We can now consider an approximate 
quantization of this system where onoy the modes of the field that carry 
infinitesimal energy are excited. This is just the quantum mechanics of n 
identical particles on the plane ( The classical thermodynamics of this sys- 
tem has been studied in |20fl .) In the simplest quantization scheme, we can 
assume these particles are bosons. There is no potential energy between 
them, if A has exactly the critical value. 

However, as we saw earlier, there might be more subtle interactions that 
arise from the boundary conditions on the wavefunctions; the boundary being 
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the region where a pair of particles come together. Such contact interactions 
seem to have been ignored in the literature on this subject:they arise if we 
approach the limit A = 1 from below rather than above. These lead to bound 
states of vortices, with an energy that is determined by a new short distance 
scale fj, exactly as discussed previously. In fact it is also possible to have a 
condensate of a large number n of vortices with a mean density of vortices 
(or magnetic flux) given by the function |"u(V)| 2 determined in the last sec- 
tion. It is of much interest to serach experimentally for such a condensation 
of vortices in supercondutors which are on the borderline between type I 
and type II. The prediction that the energy of such a configuration depends 
exponentially on the number of vortices also should be tested experimentally. 

6 Scalar Field in Three dimensions: a Cau- 
tionary Tale 

Now let us see why the above approach cannot work as it stands in three 
dimensions. Although a renormalized resolvent and a principal operator can 
be constructed, the spectrum is not bounded below for more than two parti- 
cles. (However there may be other aproaches that give a sensible formulation 



of this problem. See e.g., [18 



We get again, a Krein formula: 

R(E) = + b j } r b . (173) 

K J H -E H -E $0u, E) H -E K ' 

The eigenvalues are given by <&(E)\i[) >= 0. Again, we can rescale all the 
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momenta by y/\E\ to get 
where 



167T 



+ V\ E \ w 



(174) 



W 



1 

8^ 



[dpW(p)V 



H + u(p) + l 



x(p) 



- J [dp 1 dp 2 dp' 1 dp' 2 \x 1 (pi + P2) 



1 



H + w(pi) + u(p' 2 ) + w(pi) + w(p 2 ) + 1 
+A(2n) d 5( Pl -p' 1 )^(p' 2 ) ' 



Xbi+P 2 )- 



#0 + w(Pi) + ^(p 2 ) + ^(P2, 



<P(Pl)<P(P2) 

(175) 



This time the new eigenvalues of the interacting system are given in terms 
of the eigenvalues w n of W by 

2 1 



E„ 



/i 



(176) 



16nJ u> 2 

The operator W itself depends on no parameters: it determines the spectrum 
of the interacting theory in terms of the parameter /j,. Our system will have a 
finite ground state energy if and only if W is strictly positive. The eigenvalue 
closest to zero of W will give the ground state energy by the above formula. 

So far everything looks fine: it looks like the same approach as in two 
dimensions is going to work. But we will now see that in the three body 
sector the operator W has a zero eigenvalue which mens the energy of the 
ground state is infinite. 

In the sector of interest we have just one angel and one boson. The states 
will have the form 



\ip >= / [dp][dk]Tp(p\k)(j)\k)x\p)\0 > 



177) 
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Only the first two terms in the Principal Operator will constribute to this 
sector. The eigenvalue problem for W becomes an integral equation for ip: 

- 2 - ip(p + k-k'\k') 



^V(j+k 2 +^(p\k)-J[dk') T 



= I* 

+ k 2 + 2k' 2 - 2p ■ k' + p 2 16?r v/| E 



i>(p\k) 
(178) 



The total momentum p + k of the angel and the boson is conserved. By going 
to the center of mass frame, we impose p = —k: 



{l + h 2 ]u(k)-8n [ldk']j " (//) 



/'< 



-u{k) 



(179) 



2" J " V ^ ~" J l ~'"'± + k 2 + k> 2 + k-k> y/\E\ 

Moreover it is reasonable to expect that the ground state will have zero 



angular momentum: u(k) = u(k). Performing the angular integration gives 



An [dk'}- 



u(k') 



k 2 + k' 2 + k- k! 



An 



(27T) 



(2vr) / k' 2 u(k')dk' 



-l 

1 

7T JO 



^(cos^) — 



k 2 + k' 2 + kk' cos 6 



k'u(k') 



dk' In 



l + k 2 + k' 2 + kk' 
\ + k 2 + k 12 - kk' 



Then, putting 



we get 



v(k) = ku(k) 



l + -x z 



2 2 f°° 
v(x) / v(y) k>£ 

TV JO 



\ + x 2 + y 2 + xy 



dy 



/' 



V\E 



(180) 



■v{x). (181) 



. | + x 2 + y 2 — xy. 

We will now ask whether the ground state energy determined by this 
integral equation is finite; in other words if there are solutions to this equation 
for large \E\. We use the same ideas as the corresponding argument in two 
dimensions and see why they break down |l7j . |15[ . 



r^s 



We rewrite the equation as 

v(x) = 
where 

U E (x,y) = ■ 



U E (x,y)v(y)dy 



1 



TV 



^ + v/(l + |x 2 ) 



log 



r 1 
2 



X 



V +xy 



\ + x 2 + y 2 



xy 



182) 



(183) 



V\E\ ' V '■■ >■ 

The question is whether this equation has a normalizable solution for large 
|.E|; if such a solutions exists, the three body problem has a divergent ground 
state energy even after our renormalization. 

The idea in two dimension is to show that the Hilbert-Schmidt norm 
1 1 C/^7 j j 2 of this integral kernel is less than one as |.E| becomes large . Then 
the equation v = K E v will have no solution for large \E\. But in our case 
this will diverge. In fact, 



\K 



E\\2 



JO 



K E (x, y) 2 dxdy 



Consider first the integral over y: 



9{x) 



log 



\ + x 2 + y 2 + xy 



\ + x 2 + y 2 — xy 



dy. 



We will first show that there exists a constant C such that 



g{x) 



Cx 



for large x. By scaling y — >■ xy this equivalent to showing that 



log 



i? + l + y 2 + y 



2.T 



^ + i + y 2 



for large positive x. Now, 

^ + l + y 2 + y^ 2 



log 



K + 1 + y 2 - y 



2x 2 



dy 



log 



dy^C 



l + y 2 + y 



(184) 



185) 



186) 



187) 



i + y 2 



y 



dy = C (188) 
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The integral converges (in fact has a value of about C ~ 3.489). Thus we 
see that 

| \K E \ | 2 = / — -2 g(x)dx (189) 



7T Z 



/' 



+ ^(l + f ^ 



But this diverges logarithmically, since the integrand goes like - for large 
x. The corresponding integral in two dimensions converges which makes the 
problem well-posed there. This is why the proof breaks down here. 

In fact we can just put \E\ — oo and see that there is a normalizable 
solution for v. It can be determined by iterating the integral equation starting 
from the trial function vq{x) = (1 + §x 2 ) -1 . We spare the reader the details. 

7 Appendix: An Elementary Formula for In- 
verses 

Suppose a self-adjoint operator X : Ti — ► 7i can be split into 2x2 blocks, 

fa 6t\ 

with respect to a splitting H = Hi © 7i 2 

a:Wi-^Wi, b^:n l ^n 2 , b:H 2 ^H u d : H 2 ^ H 2 . (191) 

Let the inverse be split similarly: 

A ^i, a J- (i92) 

Then, 

a/?t + ftttf = o, ba + rf/3 = (193) 
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and 

aa + tffi = 1, bft + dS = 1. (194) 

Solving the first pair of equations we get two expressions for /?: 

_ Sba' 1 = (3 = -d~ l ba. (195) 

Putting either of these into the equation for a gives two expressions for it: 

a' 1 + a- l b ] 5ba- 1 = a =[a- tfd^b]- 1 . (196) 

Similarly, 

d- 1 + d^batfcT 1 = 5=[d- barW]- 1 (197) 

By combining the two ways of writing a and S we get 

a = cT 1 + aT l b ] [d - ba^b^ba' 1 , (198) 

S = d- 1 + d- l b[a - b ] d- l b]- l b ] d- 1 . (199) 

These identities will be used in the text. 
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